Robust gates for holonomic quantum computation 
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Non Abelian geometric phases are attracting increasing interest because of possible experimental 
application in quantum computation. We study the effects of the environment (modeled as an 
ensemble of harmonic oscillators) on a holonomic transformation and write the corresponding master 
equation. The solution is analytically and numerically investigated and the behavior of the fidelity 
analyzed: fidelity revivals are observed and an optimal finite operation time is determined at which 
the gate is most robust against noise. 
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I. INTRODUCTION 

One of the possible alternatives to quantum informa- 
tion processing 0, |^ , which is attracting increasing in- 
terest, is based on geometric interferometry 0, 3- In 
this case, the transformations needed to implement the 
quantum gates are realized by making the Hamiltonian 
of the quantum computer dependent on a set of control- 
ling parameters which describe suitable closed loops in an 
associated parameter space. In the adiabatic limit the 
dynamical contribution to the evolution can be factor- 
ized and the quantum gate only depends on the topolog- 
ical structure of the manifold. This should be contrasted 
with the dynamical approach to quantum computation, 
where the desired phase factors in the quantum gates 
are of dynamical origin. Geometric quantum computa- 
tion has been formulated using both Abelian Q and non- 
Abelian Q holonomies. In Ref. Q the first experimental 
demonstration of geometric quantum gates using nuclear 
magnetic resonance was presented. Since the appearance 
of the original proposals, several studies addressed the 
implementation with quantum optical M, superconduct- 
ing and semiconducting systems [lOj. 

As is well known, decoherence is detrimental for quan- 
tum computation. Despite the large body of knowledge 
accumulated to study decoherence in open quantum sys- 
tems (llj, the study of geometric phases in the pres- 
ence of decoherence and dissipation has started only re- 
cently, although with a few exceptions, and was certainly 
prompted by the interest in quantum computation. To- 
gether with many common features with the theory of 
open quantum systems, the analysis of decoherence in 
geometric interferometry rises several distinct issues that 
are of interest both as fundamental questions in quantum 
mechanics and in quantum computation. The adiabatic 
evolution, for example, cannot occur arbitrarily slow, as 
decoherence would destroy any interference. This implies 
that the decoherence processes should be analyzed in 
close connection with non-adiabatic corrections, a ques- 
tion which is not typically present in the non-unitary 



dynamical evolution of open systems. Moreover, the pe- 
riod of the evolution fixes a new time scale, compared to 
which the different components of the bath will act dif- 
ferently. Finally, in the non- Abelian case the coupling to 
the environment may (partially) lift some degeneracy and 
therefore modify the holonomy itself. These are only a 
few examples of questions which emerge when one wants 
to study geometric phases in open systems. 

Most of the attention on the properties of geometric 
phases in the presence of coupling to an external bath 
has focused on the Abehan |il [H El El El El El 
ITgl I20I |2]| case. There are however a few important ex- 
ceptions where non- Abelian holonomies in open systems 
have been investigated as well 0, 0, |2j, . Solinas et 
al. studied the influence of parametric noise on the 
scheme for holonomic quantum computation discussed in 
Ref. . Parodi et al. analyzed the effects of differ- 
ent spectral densities of a quantum thermal bath on the 
efficiency of this scheme. The quantum jump approach 
was applied by Fuentes-Guridi et al. "24] in order to un- 
derstand under which circumstances holonomic quantum 
computation is robust against decoherence. Very recently 
Sarandy and Lidar '25^1 analyzed non- Abelian holonomies 
for open systems, starting from an analysis of the mas- 
ter equation in Lindblad form for the reduced density 
matrix. 

The aim of this work is to study the non-adiabatic 
dynamics and the effects of quantum noise on the setup 
proposed by Duan, Cirac and ZoUer 0. We will present a 
class of 1-qubit holonomic quantum gates (which includes 
the NOT gate) that are intrinsically robust against any 
type of noise. The only requirement is that the noise be 
sufficiently small, so that a master equation can be writ- 
ten. The above mentioned robustness is a consequence 
of a peculiar property of this class of gates, namely the 
possibility to realize in the noiseless case a perfect gate 
transformation (i.e. with fidelity one) in a finite time. In 
particular this class exhibits fidelity revivals, that con- 
sist in an infinite number of (almost periodic) time val- 
ues at which the fidelity reaches unity. The first revival 
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is the optimal operational time for a nonadiabatic gate 
in presence of noise, because it represents the point with 
the highest fideUty among all revivals, not to mention 
the corresponding adiabatic gate, whose fidelity is far 
lower. In this respect our analysis consists in a gener- 
alization to non-adiabatic holonomic quantum computa- 
tion. Non-adiabatic holonomies have been discussed by 
Anandan ^26i] . generalizing to the non-Abelian case the 
work of Aharonov and Anandan Very recently, the 
use of non-adiabatic phases has been discussed in the 
framework of geometric computation as a way to further 
protect the computer from decoherence . In the same 
spirit we discuss this possibility in the non-Abelian case. 
We derive a master equation for the reduced density ma- 
trix of the system in the presence of a bath in the weak 
coupling approximation. This equation is numerically 
and analytically solved, displaying revivals of the fidelity 
and the existence of the aforc-mcntioncd "optimal" fi- 
nite operation time at which the detrimental effects of 
decoherence are minimized. 

This paper is organized as follows. In Section ^ we 
briefly review the concept of holonomy in the absence 
of the environment and introduce notation. In Section 
|n]we focus our attention on the specific physical sys- 
tem discussed in and study the role of non-adiabatic 
effects, which turn out to be important when discussing 
the realistic case of a finite time evolution in the pres- 
ence of the environment, which we present in Section lWI 
We derive here a general master equation for time depen- 
dent Hamiltonians and, after specializing it to our case, 
we numerically solve it in Section The fidelity of the 
operations in the presence of noise will be discussed in 
detail, using the ideal case as a reference. Conclusions 
and further perspectives are presented in Section IVII 



by the Schrodiger equation 



II. HOLONOMIES 

Let us introduce notation. Suppose that a system, gov- 
erned by a non degenerate Hamiltonian that depends on 
time through a set of parameters, evolves adiabatically, 
covering a closed loop in the parameter space. Berry |29l | 
discovered that at the end of the evolution the final state 
exhibits, in addition to the dynamical phase, also a geo- 
metric phase, whose structure depends only on the topo- 
logical properties of the manifold on which the system 
has evolved. 

The situation changes if the Hamiltonian possesses 
some degeneracies. In this case a loop in the parameter 
space realizes more complex geometric transformations 
[30j . Let us assume that the system eigenspaces (indexed 
by m) are degenerate and denote by \mk{t)) their set of 
instantaneous eigenstates (the degeneracy index k rang- 
ing from 1 to iV,„). The instantaneous eigenstates form 
an orthonormal basis 

{m'i,,{t)\mkit)) =5„im'5kk'- (2.1) 

The time evolution of the quantum system is governed 



i-m)) = H{t)m)). 



(2.2) 



with H depending on t through a set of parameters (t) . 
Having in mind quantum computation applications we 
will suppose that the family of Hamiltonians H{x{t)) is 
iso-degenerate, i.e. that the dimensions of its eigenspaces 
do not depend on the parameters and that its eigen- 
projections Pm{x{t)) have a smooth dependence on t 
(at least twice continuously differentiable) . In particu- 
lar, this implies that there is no level crossing between 
different eigenspaces. At each time t, H{t) can be de- 
composed by using its instantaneous eigenprojections 
Pm{t)^j:k\^k{t)}{mk{t)\: 



H{t)^Y.^r,^{t)Pm{t)- 



(2.3) 



We define the operator i?, whose action transports every 
eigenprojection from to to 

R{t, to)Pm{to) - Prn{t)R{t, to). (2.4) 

Its generator is hermitian (D{t, to) = D''{t, to)) and reads 

Dit,to) - -iRit,to)^^^R{t,to). (2.5) 

In the interaction picture defined by the operator R we 
have 

i?(i,io) =i?^(i,to)-ff(0-R(^,^o) = ^em(t)P™(to). 

m 

(2.6) 

and the evolution operator can be written as 



U{t,tQ) ^ R{t,tQ)Te^p\-i / (His^to) + Dis,to))ds 



(2.7) 

T being the chronological product. In the adiabatic limit 
the evolution of the state remains confined in the de- 
generate eigenspaces. The above evolution operator be- 
comes block-diagonal and, in the case of cyclic evolution 
{Prn{t) = Pm{to)), rcads: 

U{tM) ^ E^™(^0)^"''''''' '-^"-'^"'UZPrndo), (2.8) 
in 

where the geometric evolution 



f/a"=Pexp 



A™ (a;) 



(2.9) 



is given by a path ordered integral of the adiabatic con- 
nection A™(x) = Y.^^'^dx^', with 

^"(^(^)) = PMto))R\x{t),x{to)) 
d 

X —R{xit), xito))PMto)). (2.10) 
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The holonomy thus obtained is the fundamental ingredi- 
ent for reaUzing complex geometric transformations. In 
the following section we will give an explicit example. 

In the following we will also take into account nona- 
diabatic effects under the simplifying assumption that 
the eigenvalues Cm are time-independent and the con- 
nection l|2.5|l is piecewise constant. Then the operator 
D{t,to) ~ D{to,to) does not depend on time during the 
evolution Vs e [<, ^o]- Moreover, H{t,to) = H{to) in 
(|2.t)|) and Eq. H2.7|l reduces to the useful expression 

U{t, to) = e'^*"*")^^*"'*"^ g-j(t-to)(H(to)+r>(to,to))^ (2.11) 
which, in the adiabatic limit, becomes 

U{t,to) ^ ^P„We-^(*-*«)[^(*«)+^"(*«)lp„(io). 

m 

(2.12) 

We will see that for a large class of gates it is possible 
to evaluate exactly the time evolution, including all non- 
adiabatic effects. The analysis of the evolution operator 
will enable us to find an optimal working point where the 
gate is robust against noise. This optimal time is related 
to revivals of fidelity, i.e. (finite) values of time at which 
the fidelity goes back to 1. 

III. FREE IDEAL EVOLUTION FOR A TRIPOD 
SYSTEM 

A. Preliminaries 

We consider the system introduced in 0] for holonomic 
quantum computation: see Fig. ^ where three degener- 
ate levels are connected with a fourth one by Rabi os- 
cillations. The adiabatic evolution of this system was 
analyzed in several papers for different experimental im- 
plementations i?!, ,9„ ,^J. Here we review the ideal 
noiseless case, taking into account also non-adiabatic ef- 
fects that are important in the presence of decoherence, 
when the loop cannot be completed in an arbitrarily long 
time. At time t = the logical states and 1 are encoded 
respectively in |0) and |1), while \a) is an ancilla state 
used as "buffer" during the evolution. The Hamiltonian 

I 



\Do{t)) 

mt)) 

\D±{t)) 



The computational space (belonging to the degenerate 
eigenvalue 0) is 

Cs = Span{\Do{t)),\D,{t))}, (3.5) 



of the system reads 

H{t) = |e)(rio(t)(0| + niit){l\ + na{t){a\) + H.c. (3.1) 

where ^j(t) represent the time dependent Rabi frequen- 
cies of the transitions. The loop in the parameter space 
is obtained by varying flj{t) {j = 0, l,a). In our calcu- 
lations we consider ftj{t) e M, Vt. The eigenvalues of the 
system are 



{0, ±V^o{t)^ + ^l{t)^ + ^a{ty = ±n}, (3.2) 

where is 2-fold degenerate, corresponding to a 2- 
dimensional (computational) eigenspace, and Q is kept 
constant. Therefore, the parameter space is the 2-sphere 

e 




Ola 

FIG. 1: (Color online) Scheme of a tripod system: and 
1 are computational levels, while a is an ancilla state used 
for the intermediate steps of the transformation. The three 
degenerate levels are connected with an upper level e by time 
dependent Rabi frequencies Qj{t). We also show the noise 
introduced in our analysis, that induces additional transitions 
between and e. 



of radius Q, given by {^j € R| fi^ = Q'^}. Introduc- 
ing the parametrization 

r^i = sini? cos if, r^o = slnt? sin(p, fin = cosz?, 

(3.3) 

the eigenstates take the form 



(3.4) 

I 

while \D±{t)) are the bright eigenstates belonging to ztil. 

Applying the definition given in Section ^ the 
elements of the adiabatic connection (|2.1U|) for the 



= cos ip |0) — sin(^ |1), 
= cosd sin(/?|0) -I- cost? cos(/3|l) — sinz9|a), 
= (± |e) +sint? sin(^|0) +sin?9 cos<y9|l) + cos i9\a)) / {V2) . 
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FIG. 2; (Color online) Path in parameter space for the real- 
ization of a NOT gate. The solid angle spanned during the 
evolution is 7r/2. 



eigenspace m — are 



A^p — iuy COST?, 



(3.6) 



where = -i(|i?o(io)) (i?i(io)| - |^i(io)> (^o(to)l). The 
holonomy (|2.9I) for a closed loop on a sphere for the com- 
putational space reads 



C/ad = Pe" A^<^f = e-'^y fc ™« ^''v ^ exp {iay uj) , 



where u> is the solid angle enclosed by the loop in the 
parameter space. As an explicit example (that will be 
considered in the following), let w = 7r/2, and obtain 



Ut,/2 = cxp{iay 7t/2) = iay 



1 

-1 



(3.8) 



(3.7) 



(in the basis {\Do{to)), |Z?i(<o))}), that represents a NOT 
transformation (up to a phase for the state \Do)). 

Following the discussion in the previous Section we dis- 
cuss the non-adiabatic corrections to this system. In or- 
der to use Eq. ()2.11(l . we will consider the loop shown in 
Fig. El obtained by going from the pole to the equator, 
spanning then a 7r/2 angle and finally going back to the 
pole. The solid angle enclosed in this loop is equal to 
7r/2; in the adiabatic limit (when the product fir goes 
to infinity, r being the total time of the cyclic evolution 
and O the energy of the bright states) this path yields a 
NOT gate as in Eq. (|^ . 

The first step consists in constructing the operator D 
from Eq. (|3.4|) and the definition 12.5|l : its matrix repre- 
sentation, written in the basis {|£'i(to))}i=o,i,+,-7 is 



D{t,to) = 



I ip cosd sini5/\/2 (^sini?/\/2\ 

-(^cosi? z9/%/2 ^l^h 
—Lp sin 

V sini?/V2 -i?/V2 / 



(3.9) 



r 



B. Analytical results 

One can see from Eq. H3.9() that, as far as the rate of 
change of the angles Lp and d is constant in each seg- 
ment of the path, we can use Eq. (|2.11(l to calculate the 
evolution operator along the path shown in Fig. [3 The 
complete expression is explicitly given in Appendix ^ 

A noteworthy feature of the exact expression is that 
it is factorized in three terms. In the adiabatic limit it 
simplifies to 



Til — > + 00 



/ 






—irfl 











(3.10) 



r being the total evolution time needed for covering 
the loop in the parameter space and = a^r, with 
ai = 1. This represents a NOT gate for the de- 



generate subspace and yields (fast oscillating) dynamical 
phases for the bright states. 



C. Fidelity revivals 

In order to understand how far the evolution operator 
is from the ideal one, we use the fidelity, defined as 

F{t) = Tr{aad (r)a(r)} 

= TT{U^%am:^% U./2{nr)a{0)ul/^{nT)}, 

(3.11) 

where cr(0) is the density operator describing the initial 
state, assumed to be pure, and Cad the corresponding 
operator for the adiabatic ideal evolution. The mean 
fidelity (averaged over a set of input states uniformly 
distributed on the Bloch sphere) is plotted in Fig. as a 
function of the adiabaticity parameter fir. 

F(t) asymptotically approaches the value 1 (with some 
oscillations), as expected (adiabatic limit). Interestingly, 
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FIG. 3: Mean fidelity versus the cyclic time fir (noiseless 
case), n is the energy gap between the bright and dark states, 
r is the time needed to cover the loop shown in Fig. |5| The 
average is performed over a set of initial states uniformly dis- 
tributed on the Bloch sphere. 



the fidelity is exactly one for some finite values of time, 
T = T^, that are clearly independent of the initial state. 
In this case the NOT transformation is perfect, even 
though one is far from the adiabatic regime. We now 
discuss this curious feature that will turn out to be of 
interest in the search for an optimal operation time at 
which the computation is most robust against noise (see 
Section El . 

In order to obtain a formula for the times , consider 



the operator 



(3.12) 



appearing in Eq. H3.11|l . When the three arcs in the loop 
in Fig.[21are covered in equal times, one obtains a lengthy 
analytical expression for Q, not reproduced here. In par- 
ticular, the (1,1) element reads 



An^r^+9n- cos + (^)') 



9 7r2 -I- 4 f72 T-2 



(3.13) 



and by equating the above expression to 1 one gets 



^^^16P-l , 



2n 



k e 



(3.14) 



In turn, this yields, by direct substitution in Eq. H3.12I) . 



/ 1 

1 

e-'^^^ 

\0 



\ 







(3.15) 



The first fidelity revival is obtained for k = 1 and their 
approximate frequency reads 



(3.16) 



From this equation it is clear that the periodicity of the 
maxima in Fig.|31is only apparent (the frequency depends 
on fir). On the other hand, the second term in the square 
root is very small, even for the first peak, making //fi 
approximately constant. The additional seemingly opti- 
mal points appearing in Fig. |3 (at approximately double 
frequency) are not solution of Eq. (|3.14(l . 

These revivals can be important for experimental ap- 
plications: in principle they would enable one to obtain 
a perfect NOT transformation, without reaching the adi- 
abatic regime. It is important to notice that this result 
does not depend on the initial state of the system but it is 
a feature of the chosen path: indeed the operator Q{ilT) 
does not contain any information about the initial state. 

Finally, we emphasize that similar features (and in par- 
ticular the presence of the revivals in the non-adiabatic 
regime) hold for a large class of gates. For transforma- 
tions consisting in a loop which starts at the pole, spans 
a segment on a geodesies (the equator) and goes back to 
the pole enclosing a solid angle w = 7r/2n (n G N*) there 
is a straightforward generalization of Eq. I|3.14(l : 



{2n + l)n 
2nn 



\/ IGk^n"^ — 1. 



(3.17) 



This expression is valid provided that the loop is covered 
at a constant angular speed: 



scgmcntl — V^scgincnt2 



segments 



= const. 



(3.18) 



Reversing the orientation of the loops leads to identical 
results. These general observations can be of interest for 
experimental applications. 



IV. NOISE AND MASTER EQUATION 

A physical system is never completely isolated from its 
environment and, in order to take into account the ef- 
fects produced by the latter, one analyzes the dynamics 
in terms of a master equation. In the usual approach 
to this problem one assumes that the Hamiltonian of 
the system is time independent (see for instance 32]). 
For time dependent Hamiltonians a slightly different ap- 
proach is needed. Rigorous mathematical results were 
derived by Davies and Spohn jS^I- We summarize here 
the main conclusions, providing for completeness a phys- 
ical derivation in Appendix ^ where emphasis is put on 
the physical meaning of the analysis in the context of 
adiabaticity and holonomic quantum computation. 

We consider a general Liouville operator with a time 
dependent system Liouvillian 

C{t) = Cait) + CsB - Cs{t) ®1 + 1®Cb + Csb- (4.1) 

The evolution of density operator Q{t), describing the 
system and the environment, is governed by the von 
Neumann-Liouville equation 



e(i) = /:(i) Q{t). 



(4.2) 
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We assume that there are no initial correlations between 
system (whose density matrix is cr) and bath (whose den- 
sity matrix is as) and that the latter is in equilibrium 
(e.g. in a thermal state) 



g{0) = cr(0) (g> as, 



Cb(tb = 0. 



(4.3) 



The key hypothesis in the derivation of a master equa- 
tion is that the typical timescale of the evolution is much 
slower than the timescales characterizing the bath. The 
additional hypothesis in our case is that the timescale 
related to the rate of change of the system Hamiltonian 
is the slowest timescale of our problem, due to the adi- 
abaticity of the evolution. In other words, compared to 
the bath correlation time, the evolution of Cs is always 
"adiabatic." This is assured by the condition 



T-cA < 1, 



(4.4) 



where Tc is the correlation time of the bath and the energy 
gap, A = min |e„(i) — em(t)|, characterizes the rate of 
change of Cs- Under these conditions one gets (Appendix 

El 

&it) = [Csit)+Tit)]a{t), (4.5) 
where a{t) = Ttb {(?(i)} is the system density matrix and 

r(t) = 51 ^-(*) / duTvB{CsB cxp[— £o(i) u] 

X CsBC^Y>[CQ{t)u\(JB}Quj{t), (4.6) 
Qcj(t) being the instantaneous eigenprojections of >Cs, 

it). (4.7) 

(jj 

Equation (|4.5|) is the same master equation one would 
obtain by considering Cs{t) "frozen" at time t and eval- 
uating the decay rates and the frequency shifts at the 
instantaneous eigenfrequencies uj{t) = em(t) — e„(f) of 
the system Liouvillian. 

We now turn our attention to the physical system de- 
scribed in Sec. 11111 In terms of the total Hamiltonian, 



Cit)p = ~i[HT{t),p]. 



(4. 



For simplicity we consider an environment affecting only 
the transitions between levels |0) and |e); this is enough 
for our purposes. The total Hamiltonian is 



HT{t) = H{t) + HB + XHsB, 



(4.9) 



where A is a dimensionless scaling factor introduced for 
later convenience and representing the strength of the 



noise, and H{t) is the system Hamiltonian (|3.1() . The 
bath is an ensemble of quantum harmonic oscillators. 



Hb = y^cjfcQfc^gfc, 



(4.10) 



with ujk the frequency of the fc-th mode. The interaction 
Hamiltonian is 

i^SB =^7/c(|0)(e| + |e)(0|)® (flfct+a,,), (4.11) 



where "fk is the coupling constant between the system 
and the k-th mode of the bath. By using Eq. H3.4|l we 
can write 



|e)(0| 



cos ip 



V2 

sin if cos "d 



{\D+{t))-\DAmDoit)\ 



V2 
sin if sin ?? 



i\D+it))^\D^m{D,it)\ 

i\D+{t))~\D^mi{D+it)\ + {D^m. 

(4.12) 



In the interaction picture generated by the operator R 
defined in 12. 411 . the density operator reads 

PBit) = R^(j{t)R. (4.13) 

By taking the time derivative of Eq. 14.13(1 we obtain 

a{t) = RpR{t)R^ + R[iD{t, 0), pfl(t)]i?t. (4.14) 

By plugging Eqs. H4.12|l and (|4.14() in Eqs. 1(4. 5() and (|4.6|) . 

recalling the action H2.4|l of R{t, 0) and considering that 
the eigenvalues of the Hamiltonian are time independent 
(Lamb shifts and decay rates will be time independent 
too), we obtain the following master equation: 

PR(t) = ~i[Hs[Q),PR{t)] - i[D{t, 0), pR(t)] + \^T[pR{t)] 

(4.15) 

where 

i^ap[LafiLaP^ PR{t)] 



a,/3=0,l,± 

Eq/j 



{{Lai3L^^p,PR{t)} - 2Ll^pR{t)La(3'^^ 

(4.16) 

The Lindblad operators read 

Lo.f^ = \D^m{D piO)l (4.17) 
while fai3 = //3q, with 



p cos if r sin cos i) r sin ^ sin 

JaO 9 ; Jal n ; /a± 7 



- f 

and /oo = /n = /oi = 0. 



(4.18) 
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In the case of a thermal bath the Lamb shifts and the 
decay rates read 



F 



duj- 



, - eo 
27r^th(eQ - e/3) 



(4.19) 



where P denotes the principal value, 



1 



coth + 1 



coth ^ — 1 



(4.20) 



is the thermal spectral density, ^(w) being the bare spec- 
tral density of the noise [^(w) = for w < 0)] 



^{^) = ^lk' 5{uJk - w) 



(4.21) 



and n b (w) — 1/ (exp(/3cij) — 1) the mean number of bosons 
at frequency w. 

Notice that a non vanishing temperature entails an ef- 
fective modification of the form factors of the interac- 
tion, making them in general unbounded from below |34| . 
Moreover, note that when = ep, Eqs. (|4.19l) particu- 
larize to 



lo ^ 



^^(0+) 

f3 



(4.22) 



and the Lamb shift Aaa is temperature independent. 



V. EVOLUTION IN PRESENCE OF NOISE: 
OPTIMAL WORKING POINT 

Equation H4.15|l was numerically integrated along the 
loop in Fig. 121 when the three arcs are covered at a con- 
stant angular speed. We set F+q = F+i = Fo- = 
Fi_ = 1.10, Fo+ = F_o = ri+ = F_i = 0.8 O, F 

F__ 

0.8 n 
0.70 



-+ — 

= in, F+_ = 1.2 0, F_+ = 0.7 0, A+o = A+i = 
= Ai_ = -1.10, Ao+ = A_o = Ai+ = A_i = 

A++ = A = 10, A+_ -1.2 0, and A_+ = 

These values of the Lamb shifts and decay rates 



are physically meaningful and have been chosen, some- 
what arbitrarily, for illustrative purposes. We will con- 
sider later the realistic case of an Ohmic bath at different 
temperatures and will directly derive the values of all the 
constants from Eq. H4.19|l . 

In Fig. 0] we show the behavior of the fidelity for three 
different initial states. In each graph, from top to bot- 
tom, the dissipation constant increases from = to 
0.05. 

In the noiseless case (upper line in the three plots) the 
fidelity tends to 1 when Or 00 (adiabatic limit). This 
asymptotic value is not reached monotonically: there are 
some oscillations, with maxima at _F = 1 in the noiseless 




FIG. 4: (Color online) Fidelity F versus cyclic time Qt for 
three different initial states, Q being the energy gap between 
the bright and dark states and r the time needed to cover 
the loop shown in Fig. El (NOT gate), (a) "up" state; (b) 
"down" state; (c) completely symmetric state. In each graph 
the dissipation constant A^, defined in Eq. I|4.15|l . increases 
from top to bottom: A^ = (noiseless case), 0.005, 0.01, 0.02, 
0.03, 0.04 and 0.05. 



case. This is the case discussed in Section lllll the NOT 
transformation is perfect, even though one is far from the 
adiabatic regime, at the time values given by l|3.14|l . 

Notice that although the oscillations and the general 
behavior of F depend on the initial state, the overall 
trend is of general validity and only depends on the path 
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25 50 „ 75 100 



FIG. 5; (Color online) Mean fidelity (F) versus cyclic time 
fir. All parameters are identical to those of Fig-E] 




25 50 ^ 75 100 



FIG. 6; (Color online) Mean fidelity (F) versus cyclic time 
Qt, evaluated for an Ohmic spectral density at different tem- 
peratures T = , with — 0.01. Notice the different scale 
on the vertical axis, as compared to the previous figures. 



in the parameter space and the speed at which it is cov- 
ered. Interestingly, the position of these maxima is only 
weakly dependent on noise: the horizontal shifts of the 
peaks are almost irrelevant. (This is due to the small in- 
fluence of the Lamb shifts.) The mean fidelity (averaged 
over a set of input states uniformly distributed over the 
Bloch sphere) is displayed in Fig.O One observes all the 
interesting features discussed before. 

For all the above reasons, we can estimate the value of 
the optimal gate-operation point for experimental appli- 
cations as that point corresponding to the first peak. An 
important remark is however in order: in the non adi- 
abatic regime, the gate is no longer purely geometrical. 
Both dynamical and geometrical effects contribute to the 
transformation and cannot be easily separated, because 
in general the generators do not commute. In principle 
it would be possible to extract the geometrical contribu- 
tion, but one would not gain any additional information, 
useful for experimental purposes. 



In general, in the presence of noise, the fidelity de- 
creases as the time needed for the transformation in- 
creases. This can make it difficult to obtain a pure geo- 
metrical transformation (because of the necessary adia- 
batic condition) . This problem could be slightly reduced 
by choosing a large energy gap between the degenerate 
computational space and the bright states. However, 
it seems more convenient to exploit the presence of the 
peaks and partially neglect this physical request. As a 
matter of fact, the fidelity decrease due to the noise is 
very small in the non adiabatic regime. This makes the 
NOT gate feasible by using a fine tuning of the total op- 
eration time. The best performance is obtained for an 
optimal operation time 

r*^rU-^Vr, (5.1) 

corresponding to the first peak of the fidelity. 

On the other hand, some interesting features of the 
fidelity do depend on the initial states and are apparent 
in Fig.^ Some states are less robust than others; this is 
due to the fact that, during the evolution, the population 
transfer between the levels depends on the initial state: 
thus, the longer the population "lives" in a level which is 
subject to noise, the less efficient is the transformation. 
A critical point is obviously the total amount of noise. In 
our simulations we have considered a noise strength 
ranging from 0.5% to 5%. A realistic physical estimate 
yields a noise not exceeding 0.5%. In this regime the 
fidelity at the optimal point reaches values greater than 
0.9 for all the states considered. From this result it is 
clear that we can exploit the optimal times for realizing 
the NOT transformation with a relatively high fidelity 
even in absence of additional control. 

These general conclusions can be corroborated by con- 
sidering a particular spectral density for the environment: 
Figure ini display the behavior of the mean fidelity for an 
Ohmic spectral density 

^(w) = K w exp(— cj/wc) (5-2) 

at different temperatures. We used Eq. H4.19|l and set 
K = 1/100, = 100 and T = 1//3 = O/10,f^,5f7, lOf], 
corresponding to a thermal bath at low, intermediate, 
high and very high temperatures, respectively. Notice 
that the mean fidelity is always higher than 80% in the 
whole range of times considered (up to 100 fir). In par- 
ticular, at the optimal time, fidelity decreases only up to 
a few percent, even in the very high temperature case. 

As already emphasized at the end of Sec. IIIII these 
results can be extended to more general loops, yielding 
optimal times like in Eq. (|3.17|) . This can be of interest 
for experimental applications. 

VI. CONCLUSIONS 

We studied some aspects of holonomic quantum com- 
putation by focusing our attention on a particular phys- 
ical system, shown in Fig. ^ undergoing a loop in the 
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parameter space that, in the adiabatic Hmit, yields a 
NOT gate. We derived a general expression for the evo- 
lution operator, without taking the adiabatic limit. After 
specializing this formula to a specific loop, an exact ex- 
pression for the propagator was obtained. In order to 
gain more physical insight we considered the fidelity of 
the transformation compared to the adiabatic one. We 
found that there exist some values for the duration of 
the evolution for which the fidelity is equal to 1, even 
though one is far from the adiabatic regime. The pres- 
ence of these peaks is important for experimental appli- 
cations: if the total operation time can be fine tuned to 
the first peak, one can realize a transformation which is 
the most robust against noise. In particular, we consid- 
ered several initial conditions and analyzed the effects 
of quantum noise on the evolution of the system (if the 
rate of change of the system Hamiltonian is much smaller 
than the typical timescales of the thermal bath). As an 
example, we considered a particular noise, inducing tran- 
sition outside the computational space, and obtained the 
relevant master equation. Its numerical solution yielded 
information on the behavior of the fidelity and showed 
how important the optimal points is: actually, it enables 
one to obtain (without external additional control) high 
values of the fidelity before the system suffers much from 
the detrimental consequences of the noise. It will be im- 
portant to understand whether and how this feature can 
be of help when one tries to control the system in order 
to reduce decoherence and dissipation. 

We conclude by emphasizing that the strategy sug- 



gested in this article in order to minimize the effects of 
decoherence, being based on the determination of a (fi- 
nite, non adiabatic) optimized operation time, is some- 
what different from the other strategies suggested so far 
for suppressing decoherence. Clearly, these observations 
can be of interest for experimental applications. One 
should try and understand whether the fidelity revivals, 
and consequently optimal operation times, exist also for 
two-qubit gates. 
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APPENDIX A 

We give here the exact analytical expressions of the 
operators that describe the evolution of the system in 
the loop shown in Fig. [3 A straightforward but lengthy 
calculation yields (in the basis {|£'i(0))}i=o,i,+,-) 

U^,2{nT) - [/3(r!T3) U2{^T2) Uli^Tl) (Al) 

with (we set ri — T2 — — t jZ for simplicity) 



/I 




Vo 











7r sin a 



cos a - 



iOr sin ol 



2 a 3\/2q 

\/2 (2 a n T+3 71- cos a) 3 7r-|-2 i O r cos Q:+6a sin a 
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2a 



3 a cos a—i O r sin a 

3 Vfa 
3 a cos g-h^ ^ t sin a 

3 V2a 



- 

1 




\/2 (2 2 O T+3 TT COS a) 



3 7r+2 i n T cos a+6 a sin a 

/3 



^/2 {2iQT-3Ti cos ct) \ 



-3 7r+2 i O r cos q:+6 q sin a 
/3* 







— 3 TT — 2 1 n r cos a+6 a sin a 3 tt — 2 z O r cos q+6 q sin a 



(A4) 



where a = VQtt^ + 4JFt^ / 6 and /3 = 6tt + iiftT. 
These are the operators appearing in Eq. (|3.10|) [with 
Ti = t/3, i — 1,2,3 and the angular velocities in Eq. 
(|3.11() kept constant]. 



APPENDIX B 

We derive here the master equation H4.5|l . by focusing 
on the physical aspects of the proof, in the context of 
quantum computation with holonomic gates. 

Let us start by considering the Liouville operator (|4.1() 
acting on density matrices p. Consider a density operator 
g{t) describing the system and the environment. Its evo- 
lution is governed by the von Neumann-Liouville equa- 
tion (|4.2|l . with the initial condition H4.3|l . In the inter- 
action picture engendered by the free Liouvillian £o(i), 

$o(i,0)=Texp|^ /:o(sKs|, (Bl) 

the density matrix takes the form 

giit)^^o''{t,0)g{t) (B2) 
and Eq. (|4.2|l reads 

giit) ^ CsBit,0)giit), (B3) 

where 

£ss(i,0) = $o"'(i,0)£sB*o(i,0). (B4) 
By formally integrating (jB3p we get 



= Texp|^ Csb{ 



(s,0)dsU(0), (B5) 



whose expansion up to second order in the interaction 
Liouvillian Csb reads 



1+ / dsCsB{s,0) 
Jo 

- ds duCsB{s,0)CsB{u,0) 
Jo Jo 



g{0). 



(B6) 



r 



The reduced density operator of the system is obtained 
by tracing ljB6|l over the bath: 



ai{t) = TvB{gi{t)} 



a(0)+ / dsTrB{/:ss(s,0)e(0)} 



I ds [ duTvB{CsB{s,G)CsB{u,Q)g{0)}. 
Jo 

(B7) 



The second term in this expansion vanishes for a bath at 
equilibrium (|4.3|l . Therefore Eq. I|B7|) can be written as 



ai{t) = a(0) + 



/ ds'^s~\s,Q) 
Jo 



du K{s, u) 



$5(s,0)a(0), 
(B8) 



where 



$s(i,0) = Texp|^ /:s(s)ds|, (B9) 
K{s,u) =^ Ttb ^CsbCsb (s, u)aB | (BIO) 



and 



Csb{s,u) ^ ^o{s,u)CsB^o ^{s,u). 



(Bll) 



At this point one assumes that the typical timescale of 
the evolution is much slower than the timescales charac- 
terizing the bath and derives a master equation in the 
Markov approximation. We make a further hypothesis, 
justified by the physical nature of the process we in- 
tend to study. In our system (|2.2|I - H2.3|I there is another 
timescale, related to the rate of change of the system 
Hamiltonian: we assume that this is the slowest timescale 
of our problem, due to the adiabaticity of the evolution 
and then, compared to the bath correlation time, the evo- 
lution of Cs is always "adiabatic." This is assured by the 
condition (14.411. This condition allows us to write 



$s(s,w) = Texpjy Cs{t)dtj 

~ Texpjy" £s{s)dtj ^ exp{£s{s){s - u)} 



(B12) 
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when |s — u| < Tc- Moreover the bath correlation function 
rapidly vanishes for times larger than the correlation time 

K(s,u)~0, for \s-u\>Tc. (B13) 

Using Eqs. HB10|I and HB12|I we obtain 

K{s,u) ~ Tr_B {£sBexp[£o(s)(s - u)] 

X £s_Bexp[-£o(s)(s-u)]crB} (B14) 

and thus, in the secular approximation, forced by 
Davies's projection, Q^{s)[- ■ •]Qw(s), 

r(s) = y^Qu^is) I duK[s,u) QUs) 
Uo J 

- ^Q'^i^) / duTiB {CsBexp[~CQ{s)u] 

X. CsBe-xp[Ca{s)u]aB}Qu,{s), (B15) 

where Quj{s), given by H4.7|l . are the instantaneous eigen- 
projections of £5. Therefore, Eq. (|B8|) takes the form 



a/(0-a(0)+ / ds$5"'(s,0)r(s)$5(s,0)a(0), (B16) 
Jo 



which yields the differential equation 



cri{i) = $s"'(i,0)r(t)$s(i,0)a/(t). (B17) 



Going back to the Schrodinger picture, 



a{t) = $s(t,0)a/(t). 



(B15 



we finally get 



d{t)^[Cs{t)+T{t)]a{t). 



(B19) 



This is the master equation (|4.5() : it is the same equation 
one would obtain by considering Cs{t) "frozen" at time 
i, by applying the standard Markov approximation and 
evaluating the decay rates and the frequency shifts at 
the instantaneous eigenfrequencies uj(t) = Cmit) — ^n{t) 
of the system Liouvillian. This supports our physical 
intuition of a system Hamiltonian adiabatically changing 
in a faster environment. 
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